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摘要 : ADOR 测量 是 一 种 应 用 于 深 空 探测 中 的 差分 VLBI 技术 ， 通 过 对 探测 器 和 射电 源 的 测量 
结果 进行 差分 处 理 ， 从 而 减弱 或 消除 传播 介质 、 台 站 设备 、 时 钟 等 引起 的 测量 误差 。 目 前 我 国 
ADOR 测量 的 相关 后 处 理 中 ， 有 以 首 颗 射电 源 作为 参考 源 和 邻近 射电 源 作为 参考 源 的 两 种 不 后 
时 延 修正 方法 。 利 用 对 嫦娥 四 号 实时 任务 期 间 VLBI 测量 数据 对 这 两 种 方法 进行 研究 分 析 ， 得 到 
两 者 的 时 延 结 果 相 差 在 2 ns 以 内 ， 且 存在 随时 间 变 化 的 趋势 ， 通 过 对 基线 时 延 闭 合 差 以 及 定 轨 
残 差 进行 分 析 发 现 ， 相 较 于 首 颗 射电 源 修正 的 方法 ， 邻 近 源 修正 存在 一 定 优势 。 
X 键 词 : VLBI; ADOR; 嫦娥 四 号 ;时 延 修正 
中 图 分 类 号 : P228.6 文献 标识 码 : A 
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术 提 供 测 轨 服 务 ， 最 快 能 在 40 s 以 内 即 可 完成 4 个 射电 望远镜 同步 数据 采集 、 实 时 数据 传 
输 和 对 6 条 基线 的 数据 实时 处 理 ， 获 得 高 精度 时 延 结果 
获取 深 空 探测 器 的 高 精度 位 置 和 轨道 是 深 空 探测 活动 的 关键 。 早 期 的 深 空 探测 任务 主 
要 使 用 无 线 电 测 距 和 测速 的 方法 对 探测 器 进行 测定 轨 工 作 ， 通 过 测量 信号 在 探测 器 和 地 面 
测 站 之 间 的 传输 时 间 得 到 距离 ， 利 用 信号 的 多 普 勒 频 移 获得 探测 器 的 速度 。 测 速 、 测 距 得 到 
的 是 探测 器 视线 方向 的 速度 和 距离 信息 ， 当 测 站 与 观测 目标 间距 离 增 加 时 ， 其 在 垂直 视线 
方向 上 的 精度 缺陷 ， 使 其 测量 结果 难以 满足 一 些 关 键 弧 段 的 定 轨 需求 ， 为 此 ， 美 国 深 空 网 
(Deep Space Network, DSN) 于 20 世纪 70 年 代 开始 引入 VLBI BOR" ^ 

VLBI 技术 是 广泛 应 用 在 天 体 测 量 、 天 体 物 理 等 领域 的 测量 技术 ， 是 T 前 分 辨 率 最 高 的 
天 文 观测 技术 ”。 它 利用 2 台 以 上 相距 很 远 的 射电 望远镜 同时 观测 同一 目标 ， 接 收 其 发 射 的 
电磁 波 信号 ， 通 过 测量 信号 波 前 到 达 不 同 测 站 之 间 的 时 间 差 (时 延 ) 及 其 变化 (时 延 率 )， 可 
以 提供 观测 目标 在 切 向 方向 的 位 置 与 速度 信息 。 由 于 其 在 垂直 视线 方向 上 的 灵敏 度 优 势 ， 可 
与 测速 、 测 距 数 据 互补 得 到 更 精确 的 航天 器 位 置 和 轨道 信息 ”， 目 前 该 技术 被 美 欧 、 日 本 等 
国 的 航天 机 构 作为 主要 测量 手段 之 一 应 用 于 深 空 探测 器 的 测定 轨 工 作 中 。 在 enn 测量 过 程 
中 信和 号 先后 经 过 电离 层 、 对 流 层 ， 而 后 到 达 地 面 射电 望远镜 ， 经 过 天 线 聚 集 信 号 、 接 收 机 对 
信号 放大 、 混 频 、 滤 波 处 理 后 被 数据 采集 终端 采样 记录 。 由 于 电离 层 、 ee GAA 
影响 时 延 测 量 结果 ， 而 且 通 常 探测 器 发 射 的 信号 较 弱 、 带 宽 较 罕 ， 为 提高 测量 精度 ， 和 常 需 
要 交 蔡 观测 射电 源 和 探测 器 ， 并 使 用 射电 源 测量 结果 对 探测 器 进行 修正 ， 这 种 方法 称 为 差 
分 VLBI (delta very long baseline interferometry, 和 VLBI)。 在 深 空 探测 中 ， 天 线 接收 的 是 
探测 器 发 射 的 DOR 信和 号， 因此， 这 一 测量 方式 也 被 称 为 ADOR 型 VLBI 技术 ””， 或 双 差 
单程 测 距 (delta differential one-way ranging, ADOR ). ADOR 相 较 于 其 他 测量 方式 可 以 
提供 天 平面 内 直接 的 几何 测量 ,减少 了 每 周 跟 踪 观 测 的 弧 段 数 ， 同 时 ADOR 无 需 发 送 上 行 
信和 号。 测速 、 测 距 对 探测 器 视 向 方向 变化 敏感 ， 而 VLBI 对 垂直 视线 方向 的 变化 敏感 ， 两 
者 万 为 补充 ， 在 一 些 关 键 阶段 ， 如 地 月 转移 、 近 月 制 动 阶段 更 能 体现 其 优势 ”。 

我 国 从 剥 娥 一 号 开始 将 VLBI 技术 应 用 于 深 空 探 测 器 的 测定 轨 ， 在 嫦娥 一 号 和 娣 娥 二 号 
中 采用 8 波段 遥测 信号 开展 VLBI 测 轨 工作 ， 嫦 娥 二 号 任务 中 开展 了 X 波段 ADOR 测量 技 
术 实 验 ， 采 用 事后 处 理 模 式 进 行 了 数据 处 理 分 析 工 作 。 剥 娥 三 号 任务 中 首次 实现 了 ADOR 
数据 实时 处 理 能 力 。 嫦 娥 四 号 任务 中 ， 使 用 S 波段 和 DOR 技术 对 中 继 星 开展 了 测定 轨 工 
作 。 在 嫦娥 一 号 任务 中 ， 采 用 首 颗 射电 源 修 正方 法 ， 在 嫦娥 二 号 任务 中 采用 了 邻近 射电 源 修 
正方 法 ， 即 以 1 h 为 轮换 周期 进行 探测 器 、 射 电源 的 观测 ， 其 中 射电 源 观 测 15 min， 探 测 
器 观测 45 min。 在 娣 娥 三 号 任务 中 使 用 快速 切换 观测 模式 的 邻近 射电 源 修正 方法 ， 即 射电 
源 和 探测 器 各 5 min 交替 观测 。 在 嫦娥 四 号 探测 器 实时 任务 中 ， 同 样 采 用 快速 切换 观测 模 
式 ， 完 成 了 嫦娥 四 号 探测 器 的 测 轨 任 务 。 本 文 对 ADOR 测量 在 探 月 工程 中 的 应 用 做 了 介绍 ， 
同时 为 探究 首 颗 射 电源 修正 与 邻近 射电 源 修正 方法 之 间 的 差异 ， 分 别 采用 两 种 时 延 修正 方 
法 对 娣 娥 四 号 实时 任务 期 间 的 ADOR 测量 数据 进行 了 处 理 ， 并 对 处 理 结果 进行 对 比分 析 。 
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2  ddor 测量 技术 


2.1 ddor 测量 原理 


通常 ADOR 按照 “射电 源 一 探测 器 一 射电 源 ” 的 时 序 分 别 对 射电 源 和 探测 器 进行 观 


测 。 图 四 展示 了 ADOR 测量 中 探测 器 和 射电 源 的 相对 位 置 关系 ， 以 及 两 者 信号 传播 路 径 和 
VLBI 系 统 内 部 的 数据 流 。 探 测 器 与 射电 源 信号 到 天 线 的 传播 过 程 中 受到 电离 层 、 对 流 层 的 
影响 ， 会 产生 附加 时 延 ， 另 外 ， 天 线 系统 内 部 传输 电缆 的 长 度 受 温度 等 外 界 环 境 影 响 发 生 的 


变化 以 及 时 钟 稳定 性 都 会 对 时 延 产生 影响 。 因 此 ， 探 测 器 (SC) 和 射电 源 (QSR) 的 时 延 测量 


值 可 以 分 别 表示 为 : 


SR 
TQSR = T$ ^ + ATSR 十 2 


Tsc E | AUT bg 


inst 


式 中 ， AT 是 几何 时 延 的 模型 误差 ; 7Tinst，7Tc 


QSR QSR QSR SR 

Tolk | Tatm | Tion | T E ( 1 ) 
SC SC SC SC 

Telk | Tatm | Tion F Te (2) 


k» Tatm; Tion 分 别 为 台 站 仪器 时 延 、 基线 钟 差 、 


中 性 大 气 以 及 电离 层 时 延 ， T. 表示 测量 中 的 其 他 误差 因素 。 其 中 几何 时 延 7 是 利用 测量 目 
标 和 测 站 坐标 以 及 地 球 定向 参数 等 先 验资 料 进行 计算 所 得 ， 是 时 延 观测 量 的 主要 部 分 ”。 
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由 于 对 射电 源 和 探测 器 的 观测 是 在 相同 频段 ， 同 样 通道 进行 ， 而 且 河 外 射电 源 和 探测 器 
角 距 很 近 ， 通 常 角 距 在 10 以 内 ;对 两 者 的 交替 观测 在 数 分 钟 内 完成 ， 信 和 号 在 传播 路 径 上 的 
误差 具有 很 强 的 时 空 相关 性 ， 因 此 ， 将 强 河 外 射电 源 的 残余 时 延 与 探测 器 的 残余 时 延 进 行 差 
分 处 理 ， 可 以 消除 如 台 站 位 置 、 传 输 介质 和 仪器 等 共同 因素 引起 的 测量 误差 *””。 这 里 的 
残余 时 延 是 相对 于 相关 处 理 模型 时 延 的 观测 改正 值 。 将 射电 源 的 时 延 观测 值 (Tosr) 内 插 或 
外 推 到 探测 器 观测 时 刻 ， 并 与 探测 器 观测 值 (rsc) 做 差 : 


AT = TSC — TQsR = (T$S — qe 十 (AT3° 一 A +E , (3) 


式 中 的 ATSO — Ar9SsR 项 ， 是 探测 器 时 延 模 型 误差 与 射电 源 时 延 模型 误差 的 差 值 ， 由 于 作为 
参考 源 的 河 外 射电 源 位 置 等 信息 是 精确 已 知 的 ， 可 认为 ATP = 0; 而 探测 器 位 置 是 未 知 
的 ， 其 时 延 模型 不 够 准确 ，A75° 不 为 0， 是 需要 求解 的 量 。 为 适应 定 轨 、 定 位 的 需求 ， 在 
嫦娥 四 号 的 测量 任务 中 采用 08€ ASO 作为 观测 量 输 出 ， 并 用 于 探测 器 的 测定 轨 处 理 ”。 
e 表示 式 四 ,四 中 后 四 项 差 值 ， 由 于 两 者 观测 时 间 和 位 置 并 不 严格 一 致 ， 故 依然 存在 一 些 
未 能 消除 的 误差 ， 该 项 也 是 ADOR 测量 的 主要 误差 源 。 
2.2 ”射电 源 时 延 修 正 算法 

如 图 加 所 示 ，X 波段 ADOR 包含 五 个 信号 ， 一 个 主 载波 和 两 对 侧 音 。 假 定 探测 器 多 信 
号 残余 相位 用 $C 表示 ， 对 应 信号 频 点 处 的 射电 源 残 余 相位 为 We , 天 为 入 波段 的 五 个 
ADOR 信号 。 经 射电 源 相 位 修正 后 的 探测 器 相位 为 : 


p= 4 — a, (4) 


QSR 的 计算 公式 为 : 

PR = ped — fo)-F do , (5) 
rib JS TRU es c ERG SE rU Sog IR AE, f. 为 探测 器 信号 IR, fo 
为 射电 源 参考 频率 ，d 为 射电 源 参 考 频率 处 的 相位 。 


FIPC (8 470 MHz) 


FILI (-3.85 MHz) F1R1 (+3.85 MHz) 
PIC MD F1R2 (419.25 MHz) 


2 DOR 信号 结构 


到 射电 源 修正 后 的 探测 器 时 延 ^s SCENE BRRR 
相位 也 未 进行 大 气 、 电 离 层 残 余 相 位 改正 ， 则 扣 


对 o, 进行 多 信号 综合 (带宽 综合 ) 得 
相位 未 进行 大 气 、 电 离 层 改正 ， 探 测 器 残 
掉 介 质 效应 后 的 探测 器 时 延 为 ; 
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7" 中 还 包括 了 从 射电 源 观测 时 刻 到 探测 器 观测 时 刻 期 间 由 于 基线 钟 速 引 起 的 时 延 改正 。 假 
定 基线 钟 速 用 射电 源 测定 ， 则 : 


Pay 7QSR — 7QSR | 7QSR  ， (7) 
于 是 扣除 钟 速效 应 后 的 探测 器 时 延 为 : 
T=7 — fimxAt , (8) 


At 为 射电 源 与 探测 器 两 次 观测 之 间 的 时 间 差 。 

带宽 综合 时 延 是 通道 条 纹 相 位 随 频 率 变化 的 斜率 ， 各 通道 的 条 纹 相 位 可 以 是 n x 360°, 
其 中 为 任意 整数 ， 由 此 产生 了 模糊 度 问 题 。 在 进行 带宽 综合 前 使 用 参考 通道 射电 源 时 延 
消除 模糊 度 。 以 双 通 道 为 例 进 行 分 析 ， 残 余 相 位 pi 与 残余 时 延 + 有 如 下 关系 : 


Qi —2nfiT —2"nN; , (9) 


式 中 ，NN; 是 相位 模糊 度 ， 那 么 带宽 综合 时 延 为 ; 
002—601 Nə- Nı 
2n(f — fi) i je—R ” (10) 
其 中 ， = 即 为 带宽 综合 的 时 延 模糊 度 。 为 确定 等 模糊 度 ， 需 要 获得 一 个 初始 时 延 观 
ME To 作为 参考 ， 值 可 以 由 通道 内 具有 一 定 带宽 的 信号 ， 如 主 载波 信号 进行 条 纹 拟 合 得 到 。 
由 初始 时 延 值 再 结合 下 式 进 行 计算 ， 可 得 到 一 个 外 推 相 位 $5: 


$»— Q1 +27 (f2 — fi) To . (11) 
对 外 推 相 位 与 实际 计算 的 相位 做 差 、 取 整 可 以 得 到 相对 相位 模糊 度 : 
M-N = E e (12) 


式 中 “[]” 表 示 取 整 将 式 (I2) 代入 式 (I0). 便 可 得 到 双 通 道 带宽 综合 时 延 。 

对 于 图 四 所 示 DOR 信号 ， 首 先 使 用 主 载波 信号 的 时 延 作为 参考 时 延 ， 由 第 一 侧 音 信号 
计算 出 三 通道 带宽 综合 时 延 ， 而 后 将 该 时 延 作 为 参考 时 延 再 进行 计算 可 得 到 五 通道 带宽 综 
合 时 延 。 

基于 以 上 分 析 ， 多 信号 探测 器 时 延 计算 流程 如 下 : 

(1) 根据 式 6) 计算 计算 各 频 点 处 的 射电 源 残余 相位 ， 并 用 式 名) 对 探测 器 单 信号 作 相 
位 修正 得 到 p; 

(2) 对 主 载波 信号 进行 条 纹 拟 合 得 到 单 通道 时 延 ， 作 为 初始 参考 时 延 0; 

(3) 根据 式 (E) — (I2) 进行 带宽 综合 ， 计 算得 到 探测 器 多 信号 综合 时 延 7; 

(4) 对 7’ 进行 介质 修正 和 钟 速 修 正 ， 可 得 到 探测 器 信号 综合 时 延 7: 


T — T — (rim t Tis) (rm Tin ) akt, d) 


RB, Tak 采用 式 人 四) 计算 。 


C— 
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3 WRU rnm ADOR 测量 


3.1 ”测量 系统 

嫦娥 四 号 VLBI 观测 网 ( 见 图 团 由 上 海天 马 射 电 望 远 镜 (TM)、 北 京 密云 射电 望远镜 
(BJ)、 乌 鲁 木 齐 南山 射电 望远镜 (UR)、 昆 明 凤 凰 山 射 电 望远镜 (KM) 以 及 上 海天 文 台 VLBI 
数据 处 理 和 指挥 中 心 (VLBI 中 心 ) 组 成 。 它 们 之 间 通 过 高 速 专用 网 络 连接 ， 组 成 具备 实时 数 
据 传输 和 处 理 能 力 的 深 空 测量 网 。 各 台 站 的 观测 数据 通过 网 络 传 输 到 上 海 VLBI 中 心 ， 经 过 
相关 处 理 、 相 关 后 处 理 ， 得 到 探测 器 的 时 延 和 时 延 率 ， 随 后 利用 VLBI 测量 数据 以 及 测 距 、 
测速 数据 对 探测 器 进行 定位 定 轨 。 


密云 站 (50 m) 


云南 昆明 站 (40 m) 


图 3 ”VLBI 观测 网 分 布 图 


3.2 ”观测 模式 

从 嫦娥 四 号 探测 器 发 射 ， 到 进行 月 球 背 面 着 陆 前 ，VLBI 每 日 都 会 开展 探测 器 的 实时 观 
测 任务 ， 其 中 有 15 d 开展 中 继 星 和 探测 器 分 时 观测 。 在 开始 对 探测 器 进行 跟踪 测量 前 会 对 
强 射电 源 进 行 30 min 左右 的 观测 ， 用 于 改正 观测 装置 带 来 的 误差 和 为 探测 器 观测 时 的 带宽 
综合 提供 各 通道 的 修正 相位 ;然后 再 进行 探测 器 与 其 邻近 的 射电 源 之 间 各 5 min 的 交替 观 
W. WKAR, ADOR 测量 的 完整 周期 按照 “射电 源 一 探测 器 一 射电 源 ” 的 观测 顺序 组 
成 ， 其 中 对 目标 的 一 次 连续 观测 称 为 scan， 一 个 测量 周期 中 有 3 个 scan。 探 测 器 与 射电 源 
之 间 的 角 距 一 般 在 2* ~ 15° 之 间 ， 在 射电 源流 量 满足 测 量 精度 要 求 的 前 提 下 ， 优 先 选取 最 
靠近 探测 器 的 射电 源 进 行 观测 ; 最 后 对 强 射 电源 再 进行 30 min 左右 的 观测 ， 用 于 事后 处 理 
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消除 各 个 信号 记录 通道 间 的 设备 相位 差异 以 获得 更 精确 的 观测 结果 。 


8 探测 如 


Zr 


时 间 


4 ”嫦娥 四 号 中 射电 源 与 探测 器 的 交替 观测 流程 


3.3 ”数据 处 理 

ADOR 测量 系统 的 数据 处 理 流程 如 图 日 所 示 ， 各 测 站 使 用 8 通道 、2 比特 量化 进行 采 
样 记录 的 射电 源 信 号 和 探测 器 信号 通过 网 络 实时 传输 到 VLBI 中 心 进行 处 理 。 相 关 处 理 机 对 
探测 器 和 射电 源 信号 进行 处 理 得 到 各 人 台 站 的 自 相 关 谱 及 它们 之 间 的 互相 关 谱 ， 对 探测 器 约 
每 5 s 一 次 结果 、 对 射电 源 约 5 min 输出 一 次 结果 。 之 后 经 过 相关 后 处 理 ， 得 到 探测 器 的 时 
延 、 时 延 率 数据 ， 用 于 对 探测 器 定位 、 定 轨 。 


探测 器 信号 射电 源 信号 
相关 频谱 相关 频谱 
探测 器 数据 处 理 | | 。 射电 源 数据 处 理 
探测 器 时 延 、 射电 源 时 延 、 
x 一 MER 
MER | 介质 改正 
探测 器 时 延 、 时 延 率 修正 
改正 后 的 时 延 、 时 延 率 
探测 器 观测 值 
| 输出 结果 
|! VLBI 中 心 [时 延 、 时 延 率 ] 


5 ”数据 处 理 流程 图 


在 相关 后 处 理 中 分 别 采 用 5 s 和 300 s 的 积分 周期 对 探测 器 和 射电 源 互 相关 数据 进行 条 
纹 拟 合 、 带 宽 综合 ， 获 得 各 基线 的 探测 器 和 射电 源 的 粗 残余 时 延 、 时 延 率 。 探 测 器 处 理 


先 分 别提 取 单 通道 DOR 信号 相位 ， 并 使 用 参考 源 来 校准 单 通道 相位 ， 之 后 计算 出 带宽 综合 
的 时 延 、 时 延 率 ; 然后 对 粗 时 延 、 时 延 率 进行 介质 改正 ， 扣 除 介质 对 时 延 测 量 的 影响 ， 在 事 
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后 处 理 模式 中 使 用 水 汽 辐射 计 、GNSS 双 频 接收 机 、 气 象 仪 等 大 气 参数 测量 设备 实时 测量 的 
数据 ， 计 算出 传播 介质 在 视线 方向 上 的 时 延 来 修正 其 引起 的 误差 ”， 在 实时 模式 中 使 用 模 
型 预报 的 介质 时 延 来 对 粗 时 延 改 正 。 通 过 临近 射电 源 观测 ， 利 用 射电 源 时 延 修正 探测 器 时 
延 ， 进 一 步 前 弱 了 传播 介质 对 时 延 测量 的 影响 ， 以 及 台 站 钟 差 、 设 备 时 延 对 时 延 测 量 的 影 
响 。 在 实时 数据 处 理 模式 中 ， 使 用 在 探测 器 之 前 观测 的 射电 源 结果 外 推 到 探测 器 观测 时 刻 ， 
对 探测 器 时 延 进行 修正 ; 后 数据 处 理 模 式 中 ， 将 在 探测 器 前 后 观测 的 射电 源 结果 内 插 到 
与 探测 器 观测 时 刻 ， 来 进行 修正 。 最 后 将 改正 后 的 残余 时 延 与 模型 时 延 相 加 得 到 探测 器 的 时 
延 观测 值 。 
根据 参考 源 的 选择 ， 有 两 种 不 同 的 时 延 修正 方法 : (1) 使 用 首 颗 长 时 间 观 测 的 射电 源 作 
为 参考 源 ; (2) 使 用 交 蔡 观测 中 相位 、 时 延 、 时 延 率 精度 满足 质量 要 求 的 最 近 射 电源 作为 参 
考 源 。 早 期 的 任务 中 ， 对 探测 器 的 时 延 修正 使 用 的 是 方法 一 ， 在 嫦娥 四 号 中 采用 方法 二 对 探 
测 器 时 延 进行 修正 。 下 面 根据 嫦娥 四 号 实时 任务 期 间 ADOR 测量 资料 ， 按 照 事 后 数据 处 理 
模式 分 别 用 两 种 不 同 的 时 延 修 正方 法 进行 处 理 ， 对 两 者 时 延 数据 处 理 结果 进行 对 比分 析 。 


z 


TE 
lm 


X 


S 


4 探测 器 ADOR 测量 结果 分 析 


表 四 给 出 每 次 观测 的 任务 代码 和 观测 时 长 信息 ， 第 3 列 是 观测 探测 器 的 scan 数 ， 每 次 
持续 5 min, UR 站 位 于 四 个 测 站 最 西 面 ， 其 入 场 观测 探测 器 的 时 间 较 晚 ， 三 站 可 共 视 时 便 
开始 观测 ， 故 表 中 分 别 列 出 了 三 站 和 四 站 共 视 观测 探测 器 的 scan 数 。 


表 1 CE4 探测 器 任务 观测 情况 
观测 代码 “观测 时 长 /h ”探测 器 scan ”三 站 scan 四 站 scan 


s8c08a 8 41 14 27 
s8c09a 8.5 45 13 32 
s8c10a 5 21 11 10 
s8cl4a 6.5 27 9 29 
s8cl5a 10 51 8 43 
s8cl6a 7.5 27 8 19 
s8c19a 9 26 4 22 
s8c20a 9 24 4 20 
s8c21a 9 29 4 25 
s8c22a 9 28 2 26 
s8c24a 9 41 3 28 
s8c25a 12 67 3 64 
s8c26a 8.5 33 4 29 


ik: 第 1 列 观测 代码 由 六 位 字符 组 成 :“8” 表 示 2018 F, "c" R7 12 H, *08" 
表示 8 号 ; 第 2 列 观测 总 时 长 (小 时 ); 第 3 列 是 探测 器 观测 scan 数 ， 一 个 
scan 时 长 5 min; 第 3, 4 列 分 别 表示 三 站 、 四 站 共 视 观测 探测 器 的 scan 数 。 
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4.1 ”残余 时 延 

群 时 延 是 VLBI 的 观测 量 之 一 ， 它 们 需要 通过 对 获取 的 射电 源 或 探测 器 信号 进行 相关 处 
理 得 到 ; 而 残余 时 延 是 指 相对 于 相关 处 理 模型 时 延 的 观测 改正 值 ， 即 观测 值 与 理论 值 之 差 。 
通过 相关 后 处 理解 算出 残余 时 延 并 对 其 进行 改正 ， 将 改正 后 的 残余 时 延 加 上 模型 时 延 后 可 
得 到 VLBI 时 延 的 解 算 结果 并 用 于 探测 器 的 定 轨 定位 。 

以 2018 ^E 12 H 16 日 的 一 次 观测 (任务 代码 s8c16a) 为 例 ， 此 次 观测 持续 7.5 h XR 
测 器 观测 27 个 scan。 按 事后 处 理 模式 ， 分 别 用 两 种 时 延 修 正方 法 对 数据 进行 处 理 。 图 日 为 
各 基线 两 种 修正 方法 得 到 的 残余 时 延 结 果 ， 红 色 标 记 为 方法 一 ， 蓝 色 为 方法 二 。 由 于 在 9 时 
附近 ，KM 站 设备 出 现 波动 致使 BJ-KM, KM-UR, KM-TM 这 3 条 基线 残余 时 延 在 之 后 的 
scan (紫色 方 框 内 的 scan) 中 出 现 约 10 ns 的 跳 变 。 在 发 生 跳 变 之 后 的 scan 中 ， 由 于 方法 一 
使 用 此 前 的 射电 源 作为 校准 源 校准 探测 器 相位 ， 未 能 准确 消除 跳 变 ， 其 结果 始终 含有 该 跳 
AB. 方法 二 使 用 跳 变 发 生 后 观测 的 邻近 射电 源 作为 参考 源 ， 参 考 源 的 测量 信息 中 同样 包含 该 
设备 跳 变 影响 ， 其 差分 结果 中 消除 了 这 一 跳 变 。 


N un 
8 A 
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从 图 加 可 知 ， 两 种 方法 的 残余 时 延 会 随 观 测 时 间 的 增加 逐渐 产生 差异 ，BJ-UR, BJ- 
TM, UR-TM 三 条 基线 中 这 一 现象 更 为 明显 ， 观 测 开始 时 两 者 残余 时 延 值 几乎 相等 ， 结 束 时 
出 现 约 2 ns 的 差异 (BJ-TM 基线 )。 中 性 大 气 、 电 离 层 、EOP 误差 、 台 站 位 置 误差 等 对 测 
量 结果 的 影响 都 与 探测 器 和 射电 源 角 距 有 关 ， 方 法 一 使 用 首 个 射电 源 作为 参考 源 ， 角 距 逐 渐 
增加 ; 方法 二 使 用 邻近 源 ， 角 距 近似 不 变 ， 所 以 两 者 残余 时 延 也 会 随 角 距 (观测 时 间 ) 增加 
出 现 明 显 变化 。 另 外 可 以 看 到 ， 方 法 二 的 残余 时 延 会 在 个 别 scan 中 出 现 偏离 趋势 的 结果 ， 而 
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观测 当天 探测 器 并 没有 变 轨 等 动作 ， 残 余 时 延 不 会 因为 预报 轨道 出 现 突变 ， 这 些 差 异 仅 有 2 
~ 3 ns。 

从 图 四 可 以 看 到 ，BJ 相关 的 基线 中 都 有 个 别 scan (黑色 方 框 中 的 scan) 由 方法 二 计算 所 
得 的 时 延 测量 误差 存在 异常 。 考 虑 到 时 延误 差 与 探测 器 信号 的 相位 误差 和 参考 源 的 相位 误 
差 有 关 ， 而 方法 一 中 的 误差 没有 出 现 异 常 ， 因 此 很 可 能 是 由 于 异常 scan 前 后 BJ 站 射电 源 
言 号 信 噪 比 下 降 。 有 具体 原因 还 需要 深入 分 析 ， 这 里 暂 不 做 讨论 。 


3 S K 

ei b] t HE 

P 20- emi z - x 

"m e 方法 二 z "3 2.5 

HE vvevw pagodas Ñ á- ~ L Às B b. E 

CEE ooo e A NN CETE e ite 

T 7 8 9 10 11 12 13 p 8 9 10 11 12 13 

E at ~ 05 — ' 1 

A 3 i vade | 

W Ly bd W i ; a i : r [ 

EN ,| ona ^ ag E vius A: 

i go aia SS g E 00- TELE THET 3 

3 7 8 9 10 11 12 13 E 8 9 10 11 12 13 

a : KI 

"B 0.5- prti nina TIT LI >= 0.5- —ş | 

= E i t H ig H 

x 2 ipid u 1 

is ws al W | iip ili i 

K oo WA E ua a E oo IE UM 

E 7 8 9 10 11 12 13 = 8 9 10 11 12 13 
scan 数 scan 数 


7 2018 年 12 月 16 日 六 基线 测量 误差 


4.2 闭合 时 延 

为 了 解 时 延 测量 结果 的 内 符 精度 ， 我 们 对 ADOR 每 天 观测 时 延 的 闭合 时 延 做 了 计算 。 
此 处 采用 的 闭合 时 延 的 意义 是 以 地 心 为 参考 的 同一 波 前 时 延 观 测量 ， 即 各 基线 相对 于 地 心 
时 标的 时 延 观 测量 之 和 (如 果 不 考虑 射电 源 结构 因素 闭合 时 延 的 理论 值 为 0)。 在 分 析 闭 合 时 
延 时 ， 一 般 只 需要 计算 VLBI 网 中 任意 3 个 测 站 组 成 的 基线 的 闭合 时 延 值 即 可 ““， 当 计算 
接近 于 0 时 便 可 以 认为 整个 网 的 观测 时 延 是 闭合 的 ， 或 者 说 其 系统 差 很 小 。 
图 区 图 日 分 别 给 出 了 表 四 中 各 次 任务 分 别 使 用 两 种 时 延 修 正方 法 处 理 所 得 KM, UR, 
TM 三 站 基线 残余 时 延 闭合 差 的 均 方 根 和 标准 差 。 其 中 前 者 可 以 反映 基线 时 延 闭合 差 与 0 值 
的 偏差 程度 ， 即 该 方法 的 系统 差 ， 后 者 可 以 反映 基线 时 延 闭 合 差 的 稳定 程度 ， 即 该 方法 的 稳 
定性 。 结 合 图 加 和 图 回 ， 两 种 方法 的 闭合 差 都 在 0.5 ns 内 ， 相 较 方法 二 而 言 ， 方 法 一 的 各 次 
任务 间 的 闭合 差 变化 较 大 ， 但 一 次 观测 内 各 scan 间 闭 合 差 变化 相对 稳定 。 这 一 定 程 度 上 反 
映 了 两 种 时 延 修正 方法 的 特性 : 方法 一 使 用 固定 的 参考 源 ， 一 次 任务 内 系统 差 相 对 稳定 ; 方 
法 二 使 用 邻近 参考 源 ， 由 于 参考 源 在 观测 过 程 中 会 发 生变 化 ， 其 系统 差 也 会 相应 发 生变 化 。 
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4.3” 定 轨 后 残 差 

为 了 对 VLBI 时 延 测 量 精 度 进行 外 部 检 核 ， 利 用 两 种 时 延 修正 方法 处 理 所 得 的 时 延 测 
量 结 果 进 行 定 轨 ， 并 计算 出 定 轨 后 残 差 。 图 中 展示 了 任务 期 间 各 次 观测 结果 的 定 轨 残 差 的 
RMS, HF 16 号 KM 站 的 设备 异常 引起 的 跳 变 使 得 方法 一 的 定 轨 残 差 明显 变 大 ， 其 余 各 次 
观测 定 轨 残 差 均 在 0.3 ns 以 内 ， 且 两 者 差异 并 不 特别 明显 。 
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本 文采 用 两 种 射电 源 标 校方 法 对 CEA 号 探测 器 实时 任务 期 间 VLBI 实测 数据 进行 了 处 
理 ， 通 过 对 两 种 时 延 修正 方法 结果 进行 对 比分 析 发 现 ， 两 种 时 延 修正 方法 得 到 的 时 延 测量 值 
较为 接近 ， 两 者 差异 有 随 观 测 时 间 而 增加 的 趋势 ， 多 数 情况 下 在 ns 以 内 。 对 于 方法 一 来 
说 ， 由 于 其 固定 使 用 一 个 长 时 间 观 测 的 较 强 射电 源 作 为 参考 源 ， 在 一 次 观测 任务 内 ， 其 测量 
误差 相对 稳定 且 精 度 一 般 较 高 。 但 是 如 果 在 观测 期 间 出 现 了 设备 系统 的 跳 变 则 会 导致 跳 变 
后 的 结果 无 法 准确 修正 ， 因 此 对 设备 系统 的 稳定 性 和 可 靠 性 有 更 高 要 求 。 对 于 方法 二 ， 使 用 
邻近 源 作为 参考 源 ， 由 于 其 探测 器 的 观测 时 间 间 隔 、 角 距 都 相对 较 小 ， 能 够 有 效 地 修正 介质 
影响 以 及 设备 波动 带 来 的 误差 。 但 由 于 采用 临近 射电 源 ， 其 信号 强度 随 不 同 射电 源 有 强 弱 
变化 ， 因 此 时 延 修正 精度 亦 会 发 生变 化 ， 表 现 为 在 不 同 scan 间 有 较 大 的 波动 ， 弥 散 度 较 大 。 
从 系统 运行 的 稳定 性 、 可 靠 性 要 求 来 说 ， 是 一 种 更 有 效 的 方法 。 无 论 是 从 三 站 基线 闭合 差 ， 
还 是 从 定 轨 后 时 延 残 差 来 看 ， 方 法 二 测量 结果 都 相对 较 好 。 
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Application of Two Time Delay Correction Methods in the 
ADOR Measurement of Chang’e 4 Probe 
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Abstract: ADOR is a differential VLBI technology applied in deep space exploration. 
'Through differential processing of the measurement results of deep-space probe and radio 
source supply, the measurement errors caused by the media, station equipment and clock are 
weakened or eliminated. At present, there are two different delay correction methods in the 
post-correlation processing of ADOR measurement in China, which use the first radio source 
as the reference source and the adjacent radio source as the reference source. The VLBI data 
of Chang'e 4 real-time mission are used to analyze the two methods. The difference between 
the two results is less than 2 ns, and there is a time-varying trend; Through the analysis of 
baseline delay closure and orbit determination residual error, it is found that, compared with 


the first radio source correction method, the adjacent source correction has some advantages. 


Key words: VLBI; ADOR; Chang’e 4; delay correction 


